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Abstract 

There are different approaches for diffractive photoproduction of charmonia. Recently, 
a new approach is proposed, in which charm quarks are taken as heavy quarks and the 
nonperturbative effect related to charmonia can be handled with nonrelativistic QCD. The 
interaction between the cc pair and the initial hadron is through exchange of soft gluons. The 
exchange of soft gluons can be studied with heavy quark effective theory and an expansion 
in the inverse of charm quark mass m c can be employed. In this approach a simple formula 
for the S-matrix can be derived by neglecting higher orders in m~ l and relativistic correction 
related to charmonia. The S-matrix is related to the usual gluon distribution g(x) at small 
x. This result is different than those from other approaches. Confronting experiment the 
result is not in agreement with experimental measurement because large errors from higher 
order in m~ l and from relativistic corrections. Nevertheless the ratio of cross sections of J/tp 
and ip(2S) can be predicted more precisely than cross-sections. In this letter we show that 
the ratio predicted in this approach with an estimation of relativistic corrections is in good 
agreement with the recent measurement at HERA. 

PACS numbers: 12.38.-t, 12.39.Hg, 13.60.-r, 13.60.Le 
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Diffractive photoproduction of a charmonium has been studied extensively. Experimentally 
diffractive photoproduction of J/ip and ip(2S) has been measured with increasing precision at 
HERA|IJ, |3|, while there are different approaches for the process^, [S], |5], [7|, M. Recently HI 
has studied diffractive photoproduction of ip(2S) and obtained the ratio|| 

R = = 0.166 ± 0.007(stat.) ± 0.008(sys.) ± 0.007(57?), (1) 

which is in agreement with an earlier measurement by Hl|fJ]. The measured ratio has a weak 
dependence on the invariant mass s of the initial states. Although this ratio and the dependence 
on s can be explained theoretically, e.g., with the approaches presented in || |7|, but there are 
always certain free parameters in these approaches. In this letter we will show that the ratio and 
the weak s-dependence can be well explained with the soft gluon approach proposed in |J, in 
which relativistic correction for charmonia is taken into account and predictions are made without 
free parameters. 

In previous approaches^], ||, [j], |[ the nonperturbative properties of charmonium is described 
by light-cone wave-functions, which are defined with operators classified by twists. For diffractive 
photoproduction only a leading-twist wave-function is used. It is then questionable if contri- 
butions from higher-twist wave-functions are suppressed or not, and also for consistency of the 
twist-expansion the c-quark should be taken as a light quark and its mass m c should be neglected. 
However, taking m c = will cause some problems in these approaches and hence m c is not 
neglected. Because the leading-twist wave-function is unknown, either a parameterization is in- 
troduced or it is related to the nonrelativistic wave-function of charmoina. By solving Schrodinger 
equation with a suitable potential, the nonrelativistic wave-function, then the leading-twist wave- 
function is obtained through a Lorentz boost. It is clear that certain free parameters are introduced 
for determining the leading-twist wave-function. For the interaction between the initial hadron 
and the cc-pair one usually uses the small-size color dipole approximation, in which some modi- 
fication is also introduced, e.g., in 0. Hence the color dipole approximation also contains some 
parameters, which need to be determined from other sources. 

A new approach||, which is different than the above discussed, is used to study diffractive 
photoproduction of J/ip. In this approach the c-quark is taken as a heavy quark, hence one 
can use nonrelativistic QCD(NRQCD) to describe the nonperturbative property of charmonia |IU 



by noting the fact that the c-quark in J/ip in its rest-frame moves with a small velocity v. A 
systematic expansion in v can be employed. At leading order of v the chamonium can be taken as 
a bound state of a cc pair. The interaction between the cc pair and the initial hadron is through 
exchange of gluons. These gluons are soft in the diffractive region. Because the c-quark is taken 
as a heavy quark, the heavy quark effective theory(HQET) [TT|] can be used to study the exchange 
of soft gluons and an expansion in m" 1 for the exchange can be performed systematically. This 



approach was first used to study the decay J J if) — > e + e~ + n + tt where the two pions are soft fll^l . 
Keeping only leading orders, results for the S-matrix element are obtained and they are different 
than those in [§], ||. The reason for this is discussed in detail in ||. 
Considering the process 

7 + V -> V + J/ip (2) 
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in the diffractive region, the differential cross-section in the approach in |J is given by: 
da AQ 2 a e m 3 v^irp ,2 f-, ^z^- 



rft 9 s 2 ml 



(JMO^S^Jm ■ J2\Tr\ 2 ■ (l + + 0(t; 2 )) , (3) 



where X) is the summation over spin of the final hadron h and the spin average of the initial hadron. 
t is the squared momentum transfer between protons. The matrix element (J / ip\0\i^ S\)\J / ip) is 
defined with NRQCD fields ||10|| and is related to the leptonic decay width 

r(J/V> - e + e") = ocl m Ql^ 2 ■ (J/^\0 1 ( 3 S l )m) ■ {l + 0(v 2 )} . (4) 

The nonperturbative effect of the initial proton is represented by Tr and Tr is defined by a matrix 
element of field strength operators, which are separated in the moving direction of J/ip in the 
space-time. In the case s 3> m c , the contribution from leading-twist operator is dominant and Tr 
can be related to the usual gluon distribution g(x) by: 



T R \t^o « -in 2 m c a 2 s (m c )g(x c 

with 



1 — % tan(-a7r) 



(5) 



2m 2 

x c = (6) 

The term with tan(|a7r) is obtained by assuming xg(x) — > x~ a for x — > 0. In Eq.(5) we have set 
the renormalization scale fi to be m c . The total cross section is obtained by: 

^^) = f^| M .exp-«l.l (7) 

where the slope parameter b is measured in experiment for J/ijj and for ip(2S). 

It should be noted that at the orders we consider J/ip has the same helicity as that of the 
initial photon and the final proton also has the same helicity as that of the initial one. In Eq.(3) 
and (4) we also give the orders of theoretical errors. One of errors is from the expansion in m" 1 
for the exchange of soft gluons, which is at the order of \QCD/ m c- The parameter Xqcd is the 
typical scale of nonperturbative QCD and is at order of several hundreds MeV. Another is the 
relativistic correction at order of v 2 , because we only used leading order in the expansion in v for 
J/ip- These errors are likely very large and result in that the predicted cross-section of J/ip is not 
in agreement with experiment. For T, these errors are significantly smaller than those for J/ip, 
and the predicted cross-section of T agrees with experiment fairly well||. Although these errors 
are large for charmonia, but some of them is cancelled in the ratio defined in Eq.(l), e.g., the error 
at order of \QCD/ m c- Using the above results we obtain the prediction: 

In this ratio some corrections at higher orders of a s are also cancelled. Using the experimentally 
measured slop parameters bj^ = (4.99 ± 0.13 ± 0.39)GeV 2 and b^s) = (4.31 ±0.57 ±0.46)GeV 2 , 
and the experimental data for the leptonic decay widths we obtain: 

R w 0.47. (9) 
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Figure 1: Prediction of Eq.(12) for the ratio R, where the cures are drawn with different parame- 
terizations of the gluon distribution. Experimental data are from 

This is almost three times of the experimentally measured. With the prediction the ratio in Eq.(l) 
should have no s-dependence. Experimentally the dependence is observed and is parameterized 
by R oc s s with 5 = 0.24 ± 0.1 7||. This is also in conflict with our prediction. 

It is clear that the relativistic correction plays an important role for predicting the ratio. The 
correction is unknown and requires to be analyzed. The analysis for this correction will be in 
general complicated. Nevertheless, in the framework of NRQCD [|i~0f the differential cross-section 
can be written as: 

da{J J^ [ t _> = A± ■ (J/ij\0 1 ( 3 S 1 )\J/ij) + A 2 ■ (J/VlPi^OI JM + 0(v% (10) 
at 

where A\ can be read from Eq.(3) and A2 is unknown. The relativistic correction is characterized 
by the matrix element (J / ip\Pi( 3 Si)\ J/ip), whose definition can be found in flO| . With equation 
of motion it can be shown HIT 



(J/iP\P l ( 3 S l )\ J/1P) = m c {M m - 2m c ) ■ {J /i,\0 1 { 3 S l )\J/^j) ■ {l + 0(v 2 )} . (11) 

Replacing J/ip with ip(2S) one obtains the relation for matrix elements of ip(2S). It should be 
noted that the pole mass m c is not well determined yet, it lies in the range from 1.3GeV to 1.7GeV. 
Based on the relation in Eq.(ll) an estimation of relativistic corrections can be made if one takes 
m c is free parameter and sets 2m c = Mj/^ for J /if) and 2m c = M^s) f° r V ; (25') respectively, 
then the error due to relativistic correction at order v 2 is absent formally and the correction is 
taken into account. With this setting it has been shown that one could reduce the violation of 
the famous 14% in the case of radiative decays into ?/ of J /if) and V ; (2S')[0|. However, it should 
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be noted that this should be regarded as a phenomenological estimation, a detailed analysis and 
a precise determination of m c is needed to study the correction in a consistent way. 
Using the setting we obtain the ratio as: 



_ bj,^ T{ifj{2S) -> e+e-) a s (\M^ (2S) ) ( Mj^ 
1 1 ~~ * 



^(25) T(J/if) ^ e+e~) a s (\Mj/^) \ M <M2S), 

with 



x x g{x x ) 



12) 



_ ip(2S) _ J/-0 /-|o\ 

X2 - —2^—. ^ - ( 13 ) 
In Eq.(12) we have neglected the terms with tan(|o;7r). This term has little effect on the ratio. 



With available gluon distributions [fTS], |16| , [17] we obtain numerical results given in Fig.l, where 
the experimental data in |3[] is also drawn. 

From Fig.l we can see that the prediction based on Eq.(12) agrees with experimental data 
very well. Also, the weak s-dependence is predicted in which R increases with increasing s. This 
agreement shows clearly how important the relativistic correction for charmonia is. From Eq.(12) 
one can see that the most important factor is (Mj/^/M^s)) 5 = 0.419 for reducing the predicted 
R- value in Eq.(9). Hence, the relativistic correction in R can be larger than 50%. It should be 
noted that by taking 2m c = Mj/^ for J /if) and 2m c = M^s) f° r tfi^S) respectively the uncertainty 
from 0(\QCD/ m c) in Eq.(3) will be not cancelled in the ratio defined in Eq.(l), but it appears 
in the ratio at order of \qcd{M^ 2 s) ~ Mj/^) / (Mj/^M^S)) ~ 0.03 by taking \ QCD = 500MeV, 
hence the uncertainty is essentially smaller in the ratio than that in the differential cross-sections. 

For T-systems, if we neglect the relativistic correction and correction from higher orders of 
mZ l , we obtain the ratio 



<r(r(2S)) r(T(25) -> e + e 



<t(?(1S)) r(T(LS) -> e+e- 



0.39, (14) 



where we assume the same slope parameters for T(15) and T(25). If we take the relativistic 
correction in analogy to Eq.(12) into account, i.e., by taking 2rrib = Mx(is) ior T(IS') and 2m\ ) = 
Mr(2S) f° r T(25) respectively, we obtain numerical results for the ratio with available gluon 
distributions |15|, [16], [L7|]. The results are given in Fig. 2. From Fig. 2 we can see that the predicted 
ratio has a significant deviation from the one given in Eq.(14). Again, the most important factor 
here for reducing the value given in Eq.(14) is (M T(15) /M T(2 s)) 5 = 0.75. This indicates that 
the relativistic correction in the ratio is also significant and this significance may be beyond the 
expectation based on that the 6-quark moves inside a T-system with a velocity v 2 ~ 0.1. The 
reason for the significance in the ratio may be explained by the following: The pole mass m& 
is known more precisely than m c . A lattice determination gives m& = 5.0 ± 0.2GeV[|n|. With 
this value and experimental data of masses of T-systems, one knows that the matrix element of 
relativistic correction for T(IS') in Eq.(10) is negative, while for T(25) it is positive, based on the 
similar relation given in Eq.(ll) for the case of J /if). It results in that the relativistic correction in 
the ratio is the sum of absolute values of the both corrections. Hence, the correction is significant. 
Experimentally, the ratio is not measured. If it is measured, then our prediction given in Fig.2. 
can be tested. 
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Figure 2: The Predicted production ratio for T-systems, where the cures are drawn with different 
parameterizations of the gluon distribution. 

To summarize: In the soft gluon approach for diffractive photoproduction of 1 charmonium 
we take c-quark as a heavy quark, this allows one to use NRQCD to describe nonperturbative 
properties of the quarkonium and to use HQET to handle the exchange of soft gluons between the 
cc-pair and the initial proton. In this letter we use this approach to predict the ratio of the cross 
sections of J/ip and ip(2S). It is shown that the prediction by adding relativistic correction is in 
good agreement with the recently measured ratio. This indicates that the relativistic correction 
is very important for charmonia. Finally, it should be noted that our estimation for relativistic 
correction is a phenomenological estimation based on theoretical results, a detailed analysis and 
an precise determination of m c is needed to study the correction in a consistent way. 
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